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ABSTRACT 

We study the properties of building blocks (BBs, i.e. accreted satellites) and surviving 
satellites of present-day galaxies using the SAG semi-analytic model of galaxy forma- 
tion in the context of a concordance A Cold Dark Matter (ACDM) cosmology. We 
consider large numbers of DM halo merger trees spanning a wide range of masses 
(~ 1 X 10^" — 2.14 X IO^^A/q). We find higher metallicities for BBs with respect to 
surviving satellites, an effect produced by the same processes behind the build-up of 
the mass-metallicity relation. We prove that these metallicity differences arise from 
the higher peak height in the density fluctuation field occupied by BBs and central 
galaxies which have collapsed into a single object earlier than surviving satellites. BBs 
start to form stars earlier, during the peak of the merger activity in ACDM, and build- 
up half of their final stellar mass (measured at the moment of disruption) up to four 
times faster than surviving satellites. Surviving satellites keep increasing their stellar 
masses rather quiescently down to z ~ 1. The difference between the metallicities of 
satellites, BBs and central galaxies depends on the host DM halo mass, in a way that 
can be used as a further test for the concordance cosmology. 
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1 INTRODUCTION 

Recent times have witnessed a growing discussion regard- 
ing the formation of present-day galaxies via the complete 
disruption of sat ellites by more mas sive central objects (e.g. 
as proposed bv ,Searle fc Zinrilll978l . SZ78 hereafter), in line 
with the prediction from the Cold Dark-Matter (CDM) 
model of a bottom-up scenario of structure formation. Those 
disru pted satellites are usually refer red to as building blocks 
(BBs; llbata. Gilmore fc lrwinlll994l '). 

The observational study of the Galactic halo by SZ78 
provided hints pointing to a hierarchical galaxy formation 
scenario. They studied metallicities of the outer Milky Way 
(MW) stellar halo and inferred that it was formed from in- 
fall of small protogalactic fragments characterised by chem- 
ical evolution processes similar to the present-day dwarf 
spheroidal (dSph) MW satellites. Although a hierarchical 
scenario implies, to some degree, the infall of BBs as sug- 
gested by SZ78, it does not imply that BBs are characterised 
by formation histories (i.e. star formation rate evolution) 
similar to those of present-day satellites, as has been sug- 
gested by Geisler et al. (2007, G07 hereafter). 

Several observational facts provide indications that BBs 



and satellites should be characterised by different formation 
histories, since the galaxy population seems to have evolved 
with time as shown by (i) measurements of the evolution of 
the stellar mass function up to z « 5 (Drory et al. 2 005), 
and of the luminosity function (e.g. iFaber et al.ll2007l ). (ii) 
the decli ne of the star format ion rate towards z = (see for 
instance iHopk ins et al.l |2005| ), (iii) the important increase 
in the frequency of galaxy mer gers at high redshif t, ob- 
served as violent starbursts (e.g. iFontana et "al]|2004l ). and 
(iv) the differences in galaxy formation epochs, where mas- 
sive galaxies a ppear to be older than low-mass objects (e.g. 
Cowic, Songail a fc Cohenlll996l ). All these arguments indi- 
cate that galaxies at high redshift are likely to be different 
than their present-day counterparts and, therefore, BBs are 
also probably different from present-day satellites. This is 
also supported by several studies of our local group which 
indicate that Galactic halo stars, thought to come from rem- 
nants of previous accretion events, are characterised by dif- 
ferent abundance patterns than present-day dSph galaxies. 
Namely, extreme ly iron poor stars ([Fe/H] < —3) are not 
found in dSphs (|Helmi et al.ll2006l ). and typical 'accreted' 
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halo stars show a-elements to iron ratios ~ 0.1 — 0.3 dex 
higher than stars in dSphs (see G07). 

On the theoretical side, the galaxy population in 
the currently favoured hierarchical cosmology, obtained 
via either semi-analytic models or gas dynamic simu- 
lations, is also obser ved to evolve in broad agreement 
with observations (e.g. 



, ,, De Lucia et al.|[2006l : iLibeskind et al.l 

l2006l : iLagos. Cora fc Padillall2008l i~Tliis galaxy evolution is 
mainly driven by a changing frequency of DM halo merg- 
ers which show s a pe ak activity at redshifts z ~ 2 — 3 
l|Okamoto et al.l |2005| ). The particular c ase of the MW 
has been studied in a ACDM scenario (iRobertson et al] 
20051: iBullock fc Johnstoiil l2005l : [Pont et al.1 l2006al . I l2006bl : 



De Lucia fc Helmill2008l ) from which it arises that surviving 



satellites and BBs are characterised by different metallici- 
ties, as a result of differences in their formation time-scales. 

In this work, we investigate this issue using a consid- 
erably larger number of MW type haloes (a total of 142). 
Additionally, with the aim to obtain clues about formation 
scenarios at different mass scales, we also extend the study 
of building blocks and central/satellite galaxy populations 
to a wide range of DM halo masses. We focus on the sys- 
tematic differences in metallicity (i.e. logio(Z/ZQ)), age and 
formation time-scales between central galaxies (those hosted 
by the largest sub-halo in a DM halo), satellites (all other 
present-day galaxies) and BBs (galaxies that have already 
merged with a central galaxy). Taking into account the as- 
sumption of SZ78, we consider that, as a result of tidal ef- 
fects, a fraction of the stellar component of the BBs con- 
tribute to form the stellar halo of the central galaxy instead 
of converging to its bulge. 

We use the semi-analytic model of galaxy formation 
SAG (acronym for 'Semi- Analytic Galaxies') by Lagos, Cora 
& Padilla (2008, LCP08 hereafter), in combination with 
a ACDM A'^body simulation characterised by a periodic 
box of 60 h*^ Mpc, with a resolution of 1.001 x 10^ h"^ M© 
per DM particle. The simulation para meters are consistent 
with the resuhs from WMAP data jSanchez et al.l |2006| . 
fJm = 0.28, f^A = 0.72 and ag = 0.9; a Hubble constant 
Ho = 100/iMpc"\ with h = 0.72). As a resuh of the 
median-size of our comoving box, the maximum halo mass 
is M = 5.34 X IO^^Mq. We complement our study by us- 
ing the non-radiative A-body/SPH (Smoothed Particle Hy- 
drodynamics) simulations o f galaxy clust ers of masses M ~ 
2.14xlOi'^M0 (|Dolag et al.ll2005t ) used in lCora et al.l (|2008l ). 
extending our dynamical range and, consequently, providing 
a better sampling of the high-mass end of the galaxy pop- 
ulation. Because of resolution constraints, our analysis is 
restricted to galaxies with stellar masses Msteiiar > IO'^Mq; 
i.e. none of our conclusions refer to very low mass ob- 
jects such as the rece ntly discovered ultra faint dwarfs (e.g. 
ISimon fc Geyi2007D . In order to distinguish results intrin- 
sic to the ACDM scenario, we also use a version of the 
SAG model where disc instabilities and feedback from active 
galactic nuclei (AGN) are switched off (maintaining the pa- 
rameters of the other physical processes fixed, i.e, star forma- 
tion, supernovae (SNe) feedback). In the remainder of this 
letter, we will refer to SAG and its modified version as mod- 
els A and B, respectively. Finally, as we are interested in the 
properties of the MW and its satellites, we mimic such a pop- 
ulation by selecting z = Q galaxies with circular velocities in 
the range — 210 — 230kms~^. Such galaxies are hosted 



' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I 




11 13 13 14 



I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 l_ 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1, 







12 13 14 



12 13 14 



Figure 1. Average age and metallicity for the satellite galaxy 
population (dashed lines) and central galaxies (solid lines), as a 
function of their host DM halo mass. Results for models A and 
B are shown in the left and right panels, respectively. Errorbars 
correspond to the 20 and 80 percentiles. 

by DM haloes of masses 10^^ M© < Mdm < 4 x 10^^ M©. In 
the remainder of this letter, the 142 haloes selected this way 
will be referred to as MW-type haloes. 



2 BUILDING BLOCKS OF PRESENT-DAY 
GALAXIES 

As a first step, we analyse the way in which galaxy metal- 
licitiefl and stellar ages depend on the host DM halo mass; 
stellar ages are calculated using the stellar mass weighted 
mean. Fig. [1] shows the results for models A (left panels) 
and B (right panels), where dashed and solid lines represent, 
respectively, the average properties of satellite and central 
galaxies for each DM halo. Model A shows a shift in galaxy 
ages between central and satellite galaxies (upper left panel), 
where the latter are usually at least ~ 1 — 2 Gyr younger. The 
magnitude of this shift depends on the DM halo mass, where 
more massive haloes show larger age differences. These are 
a consequence of the AGN feedback, which reduces or even 
stops gas cooling in massive central galaxies at low redshifts, 
thus improving the agreement of galaxy luminosities and col- 
ors with observational results (see LCP08). In model B, cen- 
tral galaxies are younger than their satellites for DM halo 
masses Mdm > 10^^ Mo as a result of the absence of this 
heating source. Both models show a clear correlation be- 
tween galaxy metallicity and DM halo mass (lower panels of 
Fig-Qi with central galaxies being more chemically enriched 
than their surviving satellites. AGN feedback is not decisive 
for the build up of the relation between metallicity and halo 
mass for both central and satellite galaxies, as becomes ev- 
ident from model B. This tight relationship is a direct con- 
sequence of the hierarchical clustering process, where more 



^ The detailed implementation of metal enrichment in SAG (i.e. 
yields from core collapse SNe (SNe CC), SNe t ype la (SNe la) 
and low-intermediate mass stars) is described inBS). 
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Figure 2. Stellar mass-metallicity relation at 2 = of the full 
galaxy population for model A, model B, and model A and B 
without SNe energy feedback (symbols are indicated in the fig- 
ure key). Errorbars correspond to t he 10 an 90 pe r centil es. The 
observational relation determined bv lGallazzi et aP (|2005| ) is rep- 
resented by the shaded area. 

massive galaxies sit in initially higher density peaks and, as 
a consequence, start to form earlier on average 0, allowing 
them to acquire higher amounts of metals. Thus, a hierar- 
chical scenario is prone to produce less chemically enriched 
satellites with respect to their central galaxies, regardless 
of the AGN feedback modeling (also suggested by Tissera, 
De Rossi & Scannapieco, 2005). Note that the differences 
in metallicity between central and satellite galaxies become 
larger as the halo mass increases. In the case of MW-type 
haloes this difference is on average ~ 0.8 dex, and can be as 
high as ~ 1.3 dex. 

The different metallicities of central and satellite galax- 
ies arise from the combination of the stellar mass-metallicity 
relation (SMMR), and the lower increase of the median stel- 
lar mass of satellites compared with that of their central 
galaxy as the mass of the DM halo increases (e.g. stel- 
lar masses for satellites are 1.03 dex lower for Mnaio ~ 
3 X 1O"M0, and 1.35 dex lower for Mnaio ~ W^^Mq, on 
average). Figure [2] shows the SMMR for the full galaxy pop- 
ulation in model A (filled circles), model B (open circles), 
and model A and B without SNe energy feedback (crosses 
and open squares, respectively). As can be seen, SNe feed- 
back is essential to provi de a good agreemen t with the ob- 
served data (shaded area, iGallazzi et aLll2005h : neglecting to 
include it produces significantly higher metallicities. How- 
ever, the increasing trend of the SMMR is obtained inde- 
pendently of the details of the heating sources that affect 
the evolution of the baryonic component, including SNe en- 
ergy feedback. The SMMR of the individual populations of 
central and satellite galaxies follow closely that of the global 
population, with only minor differences at the low mass end. 

Notice that the previous metallicities correspond to the 
full stellar content of galaxies. Therefore, we cannot directly 
compare these results to measurements of metallicities of the 
MW stellar halo. As was mentioned above, we take the BBs 
at the time of the merger as a representative population of 
the stellar halo. This assumption imposes that, in order to 

Assuming a linear evolution for density fluctuations propor- 
tional to the scale-length of the Universe, a, and non-linear fluc- 
tuation for collapsed objects oc a^. 



be able to compare the properties of both populations, our 
analysis includes all simulated galaxies which have experi- 
enced at least one merger (with a BB) during their lifetime 
and have at least one surviving satellite. 

A first comparison shows that the stellar masses of BBs 
(just before the accretion) are only slightly lower than that 
of 2: = surviving satellites, and that BBs have slightly 
higher McoidGas/Afstciiar fractions (differences of a factor of 
~ 2 for MW-type haloes). Fig. [3] shows, for model A, the 
average metallicity of surviving satellites (dashed line) and 
BBs (solid line) as a function of the host DM halo mass. Im- 
portant differences between both populations arise for DM 
halo masses Mdm ^ 6 x 10^^ M©; at lower values, the metal- 
licities between satellites and BBs are rather similar, as a 
result of the smaller differences between formation epochs 
for BBs and satellites of galaxies in lower density peaks. We 
define the formation epoch as the time when a galaxy has 
acquired 10% of its final stellar mass, thus avoiding ambigu- 
ities arising from the resolution limit. The inset shows the 
difference between average metallicities of BBs and satel- 
lites, logio((ZBB)/(Zsat)), in individual systems; the filled 
circle shows MW-type haloes for which the difference in the 
average metallicities is « 0.2 dex. This measurement is con- 
sistent with results for the MW, where the abundance pat- 
terns shown by dSphs and the galactic halo are suggestive 
of hi gher metallicities for the MW BBs (e.g. iTolstov et al.l 
12004 ). In our sample, 70 per cent of the MW-type haloes 
show logio({ZBB)/{Zsat)) > 0.0, which indicates that the 
MW stellar content would be rather typical for its host halo 
mass. For model B, we find a very similar behaviour, indi- 
cating that the difference between metallicities of BBs and 
surviving satellites is an intrinsic feature of the ACDM sce- 
nario. This claim is also supported by the similarity of our re- 
sults for different dynamical friction time-scales; we use the 
two alternative recipes implemented in SAG, as described 
in LCP08. Our results confirm tho se found for MW-typ e 
haloes in numerical simulations (e.g. [Robertson et al.ll2005l ). 
but this is the first time such an analysis is performed us- 
ing an unbiased sample of merger trees (with at least one 
galaxy merger and one surviving satellite) extracted from a 
fully non-linear cosmological numerical simulation. 

It is important to remark that our model predicts that 
the majority of the stars in galaxies hosted by MW-type 
haloes are formed in-situ since BBs only contribute with 
~ 15 per cent of the z = stellar mass. In contrast, more 
than 50 per cent of the z — stellar mass of central galaxies 
hosted by DM haloes with masses Mdm > lO^^Mo comes 
from the accretion of BBs. 

By analysing the details of their formation (i.e. forma- 
tion epoch and history), we can explain the gap in metal- 
licity between BBs and surviving satellite galaxies. We esti- 
mate their formation time-scales by taking into account the 
look-back time (LBT) to the moment in which they acquire 
10 and 50 per cent of their final stellar mass. These are rep- 
resented for model A in the left panel of Fig.|4]by dashed and 
solid lines, respectively, for BBs (blue) and surviving satel- 
lite galaxies (red), as a function of the host DM halo mass. 
As can be seen, regardless of environment, BBs assemble 
their stellar mass in considerably shorter times (~ 0.25 Gyr 
to increase their stellar mass from a 10 to a 50 per cent of 
its final value) than surviving satellites (~ 1.5 Gyr), imply- 
ing that the former will naturally show enhanced a element 
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Figure 3. Average metallicity for the present day satellite galaxy 
population in individual DM haloes (dashed line) and for the 
BBs of central galaxies (solid line), as a function of the host DM 
halo mass, for model A. Errorbars correspond to the 20 and 80 
percentiles. The inset shows the difference between the average 
metallicities of BBs and satellites as a function of DM halo mass. 
The filled circle represents MW-type haloes. 




Figure 4. Left panel: Average LBT to the formation of a 10 
and 50 per cent of the final stellar mass of BBs and surviving 
satellites as a function of DM halo mass (line types are indicated 
in the figure key), for model A. Right panel: Central halo mass 
at the redshift of merger with the parent haloes of BBs and the 
average mass of parent haloes of surviving satellite galaxies at the 
same redshift. The average redshift for the merger is shown in the 
inset. Errorbars show the uncertainty on the mean. 



abundances in comparison to the latter. During short and 
intense star formation (SF) events, such as those character- 
istic of BBs, a stellar population achieves large [a/Fe] ratios 
since stars are formed from gas polluted by SNe CC prod- 
ucts; the onset of SNe la occurs ~ 1 Gyr after the first SF 
episode, when the bulk of SF activity has already occurred. 
A detailed study of a element abundances will be tackled in 
a forthcoming paper. Additionally, BBs also finish acquiring 
50 per cent of their final stellar mass substantially earlier 
(~ 2 Gyr) than surviving satellites as a result of their dif- 
ferent initial environments. BBs are aggregated into more 
massive objects earlier than surviving satellites; i.e. their 
initial location could be part of a higher density fluctuation 
than that of satellites and, therefore, collapse earlier, a s al- 
ready discussed in the literature (e.g. iHelmi et al.]|2006l ). In 
model B, we consistently find that BBs form both earlier 
and faster than surviving satellites regardless of their host 
DM halo mass, indicating that these two effects are natural 
in a hierarchical universe. 

These results indicate that the formation epoch of BBs 
in MW-type haloes corresponds to 2 ~ 2.7 (with ~ 8 per 



cent dispersion in LBT in our model) and that they acquire 
50 per cent of their stellar content by z ~ 2.25, whereas 
MW satellites show formation epochs at z ~ 2.1 (with ~ 5 
per cent dispersion in LBT for MW-type haloes) and ac- 
quire 50 per cent of their stellar mass by z ~ 1.25. In the 
concordance cosmology, the peak in the merger activity of 
galaxies occurs at 2 ~ 2 — 3 (e.g. LGP08), which coincides 
with the formation epoch of BBs, and only partially over- 
laps the epoch when surviving satellites start acquiring their 
stellar mass. Therefore, the cosmological peak merger activ- 
ity, which is directly linked to the maximum SF activity 
in the Universe, may be partly responsible for the different 
properties of BBs and satellites and would explain the fact 
that their properties are more similar at lower host DM halo 
masses (see Fig. [J]), for which most of the SF activity occurs 
at later times when the cosmic SF activity is already in its 
declining phase. Note that at fixed DM halo mass, the pop- 
ulations of BBs on the one hand, and of surviving satellites 
on the other, show very homogeneous properties; the rela- 
tive 1-a dispersion in ages and metallicities is of the order 
of 10 per cent, or smaller. 

The right panel of Fig. |4] shows in a solid line the mass 
of the central DM halo at the moment of the accretion of 
the parent DM haloes of BBs (when the collapse of the den- 
sity peak has just occurred and before dynamical friction 
starts the merging process between the BBs and the central 
galaxy) as a function of the 2 = host halo mass. The aver- 
age DM halo mass hosting the 2 = surviving satellites at 
the same LBT is shown by a dashed line. Additionally, the 
average redshift of acquisition of the host haloes of BBs is 
shown in the inset. We can see that, at the accretion epoch of 
BBs and regardless of environment, the mass of haloes host- 
ing BBs is larger than that of the surviving satellites. This 
result proves that BBs are embedded in higher initial fluc- 
tuation peaks than surviving satellites. This way, a cosmic 
queuing picture emerges where the progenitors of the most 
massive present-day objects start their formation process 
earlier than lower mass systems. Lower present-day mass 
objects need to w ait long er to reach this collapse threshold 
(see, for instance. IPress fc Schcchtcr 1974). 

In summary, the ACDM cosmology predicts that cen- 
tral and satellite galaxies represent two distinct populations 
with different age and metallicity distributions. Even though 
satellite galaxies in general cannot be expected to resemble 
the original galactic building blocks, they can still be used 
to probe the assumptions behind galaxy formation models 
by comparing their metal content to that of the halo stars 
of their central galaxy. 

2.1 Why are all the MW satellites so similar? 

The observational result that MW satellites are similar be- 
tween themselves in terms of metallicity properties (G07), 
and different from their host - the MW galaxy - has been a 
subject of intense debate. In this letter, we suggest the latter 
is a consequence of the same processes behind the build-up 
of the SMMR; regarding the former, we find that the disper- 
sion in the stellar mass and age for satellites seldom exceeds 
23 and 10 per cent, respectively. In particular, MW-type 
haloes show dispersions below 5 per cent for the formation 
LBT and 20 per cent for the stellar mass of surviving satel- 
lites. Such narrow spreads in masses and ages also results 
in narrow metallicity dispersions within DM haloes, provid- 
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ing a plausible explanation for the similarities between the 
different MW dSphs. 



BBs to imprint their signature in the continuously evolving 
relation between the properties of surviving satellites and 
their host galaxy stellar halo. 



3 SUMMARY AND CONCLUSIONS 

In this work, we have used the semi-analytic model described 
by LCP08 in combination with a cosmological A'^body sim- 
ulation of the concordance ACDM cosmology and hydrody- 
namical simulations of galaxy clusters, to study the popu- 
lations of surviving satellites and building blocks of central 
galaxies hosted by DM haloes of various masses (model A) . 
In order to isolate the effects arising from the hierarchical 
growth of structures alone, we also considered a model with- 
out AGN feedback and disc instabilities (model B). 
Our main results are summarised as follows. 

(i) On average, massive DM haloes host more chemically 
enriched galaxies than low-mass haloes, regardless of the in- 
clusion of AGN or SNe feedback. In model A central galaxies 
are older than the surviving satellite population for all the 
halo masses considered. In model B, central galaxies of the 
most massive haloes can be younger than their satellites. 
The latter effect arises as a result of the lack of the AGN 
feedback, which suppresses the star formation in massive 
galaxies at low redshift. 

(ii) Halo stars (analysed via the stellar content of BBs) and 
present-day satellite galeixies show clear differences in age 
and metallicity; the former are found to be more chemically 
enriched than the latter, regardless of the inclusion of en- 
ergy feedback sources as AGN and SNe. We find that these 
results arise from the earlier formation epochs and faster for- 
mation time-scales experienced by BBs with respect to satel- 
lites. The epoch of assembly of BBs for MW-type galaxies 
roughly corresponds to the peak merger activity in a ACDM 
universe. 

(iii) With respect to MW-typc DM haloes in SAG, wo find 
that, in 70 per cent of the cases, the metallicities of BBs are 
higher than those of surviving satellites, as is thought to be 
the case of the MW (e.g. G07). This indicates that the MW 
stellar content would be rather typical for its host DM halo 
mass. This particular result does not depend on the imple- 
mentation of AGN feedback and disc instabilities. However, 
it should be borne in mind that the observational data may 
still be subject to selection biases as has been proposed by 
Bullock & Johnston (2005). 

Our conclusions do not depend on the detailed phys- 
ical assumptions adopted in the semi-analytic model but 
are rather a natural consequence of the hierarchical growth 
of structures. To conclude, wo remark that the differences 
between BBs and surviving satellites can be thought of as 
a consequence of the existence of a SMMR, which is found 
both in observations and in models of galaxy formation. This 
relation indicates that objects with higher stellar masses, 
corresponding to higher peaks in the density field, are char- 
acterised by higher metallicities. Therefore, it is not surpris- 
ing that central galaxies and BBs (which collapsed early into 
massive objects and populate high peaks) are characterised 
by higher metallicities than satellite galaxies. The latter only 
recently began their merging process with their host galaxy 
and come from lower peaks in the density fluctuation field. 
This explains their lower metal abundances. In this picture, 
satellite galaxies are simply waiting their turn to become 
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